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^: ..... 

^ - . . . 

^ . It is anticipated that a number of techniques to measure the top quark mass at the LHC will yield 

"T'top with uncertainties of about 0.5— 1 percent. These uncertainties are mostly theoretical; they 
are usually estimated using parton shower Monte Carlo programs whose reliability at this level of 
precision is difficult to assess. The goal of this paper is to contrast those estimates with the results 
(~| . of NLO QCD computations for a few observables, often discussed in the context of high-precision 

Qh' top quark mass measurements at the LHC. In particular, we study the NLO QCD corrections to 

^ • the invariant mass distribution of a charged lepton and a _B-meson in lepton+jets channels. In the 

' dilepton channel we investigate the invariant mass distribution of a charged lepton and a 6-jet, 

. the average energy of the two leptons and the average energy of the 6-jets from top decays. 

' 1. Introduction 

o 

Measurements of the top quark mass with highest possible precision are useful for constraining 
physics beyond the Standard Model through precision electroweak tests. In spite of spectacular 
Tevatron results P, Q, the motivation to continue top quark measurements will remain strong 
even at the era of the LHC since one would like to see a consistency between direct and indirect 
evidences for New Physics. Absence of such a consistency will be a strong indication that our 
understanding of emergent New Physics is incomplete. 

ATLAS and CMS plan to employ a variety of methods for measuring the top quark mass 
[1, 0, @, S 0, Q • These methods can be divided into two classes. The first class includes so-called 
' matrix element methods. The idea is to fit the top quark mass by adjusting its value, to best 

^ I describe various kinematic features of ti events, using squared leading order matrix elements as 

probability density functions. Such methods typically lead to very small uncertainties in the top 
quark mass because nearly all the information about the events is utilized. However, the drawback 
of these methods is that it is hard to estimate the theoretical uncertainty in the top quark mass 
obtained in that way. 

The second class includes determination of the top quark mass from kinematic distributions 
that are sensitive to the value of mt- It is interesting that, up to now, all analyses of such 
distributions including their sensitivity to the top quark mass and the theoretical uncertainty in 
mt have been performed 0, 0, 0,SS@| using parton shower event generators, such as HERWIG 
and PYTHIA0 Kdi- As the result of such studies, it is often claimed that high precision in 
the top quark mass measurement can be achieved. Unfortunately, it is not clear how reliable 
such conclusions are since, by construction, parton showers can not guarantee high precision for a 
generic observable due to their approximate nature. It is therefore important to find alternative 
ways to estimate theoretical uncertainties in the description of relevant kinematic distributions, 



o 

\6 
o 
o 



^Even determination of the top quark mass from the total cross-section for top quark pair production requires 
computation of the acceptance since the total cross-section is never measured. Such acceptances are routinely 
computed using parton shower event generators. 
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since those uncertainties impact directly the precision of mt which can be expected from those 
measurements. 

In this regard, we point out that some kinematic distributions that are expected to be used 
for the top quark mass measurements, can be computed in perturbative QCD. Typically, those 
distributions involve top quark decay products; their computation at leading order in perturbative 
QCD is straightforward. On the other hand, precision requirements on mt make it necessary 
to go to higher orders in the perturbative expansion and, in spite of the fact that NLO QCD 



corrections to tt pair production have been known for about twenty years ll|, |12| , only recently 
NLO QCD corrections to top quark pair production and decay with all the spin correlations 
included became available [ISl [ij, |l5 1 . The availability of NLO QCD corrections to realistic final 



states and observables is a necessary pre-requisite for the high-precision analyses of kinematic 
distributions relevant for the top quark mass determination. Higher-order corrections are also 
important from a more theoretical viewpoint since they allow us to distinguish between the mass 
parameters defined in different renormalization schemes. In the context of tt production at the 
Tevatron, this issue was recently discussed in Ref. [l6j . Throughout this paper, we employ the 
pole mass of the top quark, for convenience. 

The goal of this paper is to use the computation reported in Ref. as a starting point to study 
some observables relevant for the top quark mass determination at leading and next-to-leading 
order in QCD perturbation theory and to investigate their sensitivity to input parameters. We 
will study four observables in this paper. In Section [5] we discuss the invariant mass distribution 
of a lepton and a i3-meson from the top quark decay; this is a simplified version of an observable 



suggested in Ref. 17| for the top quark mass measurement. Analysis of rriBi distribution with 
high precision requires computation of NLO QCD radiative corrections to the exclusive decay 
t—^-B + X + l'^ + v which we report in Section [21 In Section [3] we investigate the invariant mass 
distribution of the fo-jet and the charged lepton, and the distributions of the sum of energies of 
the two leptons and the sum of energies of the two 6-jets in the dilepton channel. Calculation of 
those distributions is performed using results reported in Ref. We conclude in Section H) 

2. Measurement of mt in top decays to the final state with identified B-meson 



It was pointed out in Ref. [17[ that the top quark mass can be accurately measured by studying 
top quark decays to an exclusive hadronic state. For example, one may consider the process 
pp ^ (t ^ + b + J/V') + (*—>■ W~ + b) and require that the W~ decays hadronically, 

decays leptonically, J/^ decays into a pair of leptons and b decays into a lepton (inside 
the jet) as well. Then, one may use the invariant mass distribution of a J/ip and an isolated 
lepton to determine the top quark mass. The requirement of a large number of leptonic decays 
reduces the rate significantly. However, it also reduces the combinatorial background from the 
incorrect pairing of the J/ip and a lepton. In addition, since no jet measurements are involved, 
the measurement is insensitive to jet energy scale uncertainties. As the result, a very accurate 
reconstruction of the invariant mass mj/^i and the measurement of the top quark mass become 



possible. It is expected [17|, [iSl ll9| that 0(1 GeV) error on the top quark mass can be achieved 
in such measurement. The only (serious) drawback of this method is that very large luminosity 
- about 100 fb~^ ~ is required since the rate is suppressed due to all (semi)leptonic branching 
fractions involved. However, it was pointed out in Ref. [20| that the situation can be improved 
by giving up the requirement of the leptonic decay of the 6-quark. In this case, combinatorial 
background increases but remains manageable, and the luminosity needed to reach uncertainty of 
about 1.5 GeV in the top quark mass is reduced to about 20 fb~^ [20| . 

The small uncertainty in the top quark mass that can, potentially, be achieved in those measure- 
ments is very attractive. It also sets the bar for other methods of the top quark mass measurement. 



planned at the LHC. It is therefore peculiar that the analyses in Refs. [Ill, |18|, ll9| are performed 
by using parton shower event generators to describe production of top quarks and their decays. 
The uncertainty in the extracted value of mt is estimated in those references by comparing results 
obtained using different parton shower event generators, such as PYTHIA and HERWIG |JJ|] or 
even diff'erent versions of HERWIG [S] . It is possible that parton shower event generators give 
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reasonable description of the required mass distribution and that the resulting error estimates of 
the top quark mass are trustworthy. However, it is important to check this, given the potential 
importance of the top quark mass measurement. To this end, it is useful to look for alternative 
ways to describe the process pp {t ^ + b ^ + J / t/j) + {t —i' +b), to ensure that 
the current understanding of top quark decays and b ^ B fragmentation is consistent with a very 
small error on nit that is claimed to be achievable through mj/^i measurement. 

One way to achieve that is to avoid using parton showers and, instead, to compute the nij/^i 
invariant mass distribution in the process pp— {t + b —?' J/ip) + {t +b) in fixed-order 

perturbative QCD. We describe how this can be done in this Section. To simplify the problem, we 



follow Refs. [18|, [19| in that we do not include the decay of a i?-meson to a J/V' meson, but only 
consider a &-quark fragmentation into a _B-meson. This is a reasonable first step because decays 
of i3-mesons to J/V'-mesons are well-studied at _B-factories. The energy spectrum of i3-mesons in 



top decays can be computed using the 6 — > _B fragmentation function formalism 21| which allows 
systematic inclusion of higher-order QCD effects. The observable that we study in this Section is 
the invariant mass distribution of the _B-meson and the lepton from the associated VF-decay. 
The NLO QCD calculation of the _B-meson energy spectrum in top quark decays was performed 



in Ref. 22| within the & — )■ i? fragmentation function formalism. However, the results of that 
reference can not be used directly for our purpose since leptons from W^-decays were integrated 
over. Because the primary object of our study is the invariant mass of a lepton from the W decay 
and a i?-meson from the ^-fragmentation and since we would like to be able to impose kinematic 
constraints on top quark decay products, we require a calculation of the NLO QCD corrections 
that is exclusive inasmuch as the top quark decay products are concerned. 

To perform such a calculation, we employ the dipole subtraction formalism of Ref. [2^. We 
point out, however, that we have an identified hadron in the final state. Hence, care is required 



when the dipole subtraction formalism is applied. In principle, Ref. |23[ does describe the construc- 
tion of the subtraction terms for such a situation, but since we deal here with decay kinematics 
and since massive particles are involved, we adopt a slightly different approach. As our starting 



point, we take subtraction terms constructed specifically for top quark decays in Ref. [2J]. We 
modify the subtraction procedure slightly, to allow for the identified hadron in the final state, and 
obtain a fully differential description of the decay t — > Z+z^ + B + X through NLO QCD. 

2.1. Calculation of radiative corrections to t ^ l^v -\- B + X . 

In this Section, computation of the radiative corrections to the decay t — (W^ 1^1^) + B + X 
is described. We assume that the i?-meson is produced by the fragmentation of the massless b- 
quark. The W boson is on the mass shell. We denote by x the fraction of energy carried away by 
the i?-meson in the top quark rest fram^ 



2tpi 



The differential decay rate reads 



2(1 -r2)' 



^2 ^ "'w 



(1) 



dPs _ / de dr, fx\ 

X 

where ^ = 2tb/{mf{l — r^)), dF^/d^ is the differential decay rate for the partonic decay t — > 
W + b + X and D{x) is the fragmentation function for b ^ B. One can restore the dependence on 
other partonic variables in Eq. ^ because coUinear fragmentation does not affect them. Therefore, 
Eq. (I2|) is a starting point for the computation of various kinematic distributions for top quark 
decays to final states with jets and an identified _B-meson. 



^Except for the S-meson, we denote particles and their momenta by the same label. We hope that this fact 
does not cause any confusion. 
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Our goal is to compute these distributions through next-to-leading order in perturbative QCD. 
To this end, the partonic decay width dF is expanded in series of the strong coupling constant as 



dn = dr[") + dr[^^ + drf ^ + o{al), (3) 

where the three terms refer to leading order decay rate and virtual and real contributions to the 
NLO decay rate, respectively. Because dvf'^ and dF^^'' have two-body final stated, the fe-quark 
in that decay has maximal energy. This implies that dF^,"'^'' are proportional to a delta-function 
of? 

and we can write 

dx [dT^°^ +dr^P} D{x)Fj,2{{p}) 



dFl°)^dFl^)^ 5(1-0, (4) 



d.dFf)!%l!).(^^i|p^)..3(M). (5) 



Note that we introduced the "measurement function" Fj^n to indicate external constraints that 
are applied to a n-particle final state. The measurement function depends on the momenta of 
final state particles, including the momentum of the B meson. It is assumed to satisfy the usual 
requirements of infra-red and collinear safety. 

Since virtual and real corrections are separately infra-red and collinear divergent and since 
the measurement function is arbitrary, we need to set up a calculation where all divergences in 
dr(i^) and in dr(V) 

are regulated separately. We construct the necessary subtraction term below 
following Ref. [3] closely. We begin by considering the matrix element that describes the real 
emission process t W + b + g 

drf ) cx IX3I' d$(3), (6) 

where d$^'^-* is the phase-space element for W, b and g, and introduce variables z, y to parametrize 
the scalar products that involve the gluon momentum 

bg='^{l-rry; tg ^ "^{l - r^){l ~ z). (7) 

The soft limit g — > requires z — > 1, y — !■ 0. The collinear limit g 7^ 0, gb ^ corresponds to 
y — > 0. We express the fraction of maximal energy (in the top quark rest frame) carried by the b 
quark in t^b + g + W decay through z and y variables 

Using Eq. (|8]), we find that f{z, y) — 1 and /(z, y) — )■ z in the soft and collinear limits, respectively. 

The soft and collinear limits motivate the construction of subtraction counter-terms for real 
emission corrections. Taking the difference of the real emission correction and the subtraction 
counter-term, we obtain an integrable expression 

which explicitly involves the fragmentation function. We emphasize that, as with any subtraction 
method, the counter-term is evaluated for values of momenta that differ from the momenta used 



^We count decay products of a boson as a single particle. 
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in the evaluation of the matrix element. In particular, the S-meson momenta are related to a 
particular 6-quark momenta in the following way 



PB = 77 r, PB = xb. (10) 

In the soft limit f{z, y) — J- 1, z — > 1 and b ^ b, whereas in the coUinear limit, /(z, y) z, b ^ zt 
and b ^ t. It follows from Eq. ^T0\\ that pB and pB coincide in both limits; of course, this is 
an important condition for the proper work of the subtraction counter-term. For the subtraction 
matrix element |A^3p, we employ [241] 

|>t3|' = |^2pDip(2,y), (11) 

where 

Dip(z, y) = C^gl^?' (^1 - 1 - z - r;e(l - z)^ - ^) , (12) 

and 7W2 is the matrix element ioi t ^ b + W . The term proportional to 77 distinguishes the 
t'Hooft-Veltman (77 ~ 1) and the four-dimensional helicity (77 — 0) regularization schemes. We 
need to employ the t'Hooft-Veltman scheme in our calculation since this is the scheme (combined 
with the MS subtraction) in which fragmentation functions are extracted from the e^e~ data in 
Ref. [ii]. 

The required momentum mapping is constructed in Ref. 2j|; we summarize it here for com- 
pleteness. We need to map a three particle final state {t ^ W -\-b + g) onto two particle final state 
(i — )■ W + b). We require that the top momentum does not change, i = t, so that 

t^W + b. (13) 

Since W"^ — — rn^ , 14^ is a valid candidate to be the four-momentum of the W boson after 
mapping but it has wrong energy for the two-body decay. To correct for that, we can make a 
Lorentz transformation 

= Ki^wr (14) 

The matrix A''" reads [24 1 

K'^" = q'"' + smh(a:) ^^^^ _ 

+ nrTr^2^''^ 2^ 2 {tWif^W + W^t^) - m^^tn^ - m^^W'W) , (16) 
(tW)'' — m^mf 

smh{x) = ^ ^ (-{ml - ml,)tW + (m^ + ml,) J {tWY ~ ml^mf] . (17) 
Applying the Lorentz transformation to VF, we obtain a simple expression 

W^a(w-^-^t)+pt, (18) 



where 



where 



^ y/imj- raj^f - 4m^7^ (m^ - m^) 

Explicit knowledge of the matrix K^'^ is required to account for momenta changes of the decay 
products of the W boson. 

We are now in position to discuss the integration of the subtraction term over the unresolved 
phase-space. We use the phase-space factorization, described in Ref. [24], d^^'^) = d$(^Mg, where 



■ (l-r)^mr" (4^)- 

167r2 r(l - 



2^ i ymax 

(l^) /d^(^' + ^(i-^'))"' / ^/"'(y-ax-2/)-^ (20) 
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In Eq. (|20p . e is the dimensional regularization parameter and ^max ~ (l+r)^z(l — z)/(z+7'^(l — z)). 
We need to calculate 

(21) 



/dip y" ^5 Dip(z, y) z ^ ^ 



The two integrals over y that we need (through an appropriate order in e) are 

ymax , „ , Ifmax 



— 2/"'(2/max - y) 
2/ 



■y he 
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Upon integrating over y, we obtain 



/dip — 



(47r)2-2^r(l-e)(l-r2)2^ ] z \z 



'^d{^-^ {t" + z{\-T')rz-^\\-z)- 
2z{\^2e) 







1 - z 



1 - z 



1 -K' 
£ 6 



(z + r2(l-z))(l-z) 



(22) 



Because the z-dependent fragmentation function D is present in Eq. (1221) , we can not integrate 
over z analytically. The best we can do is to extract infra-red and collinear divergences. We find 



27rr(l-e) V "1? 
I>(.)(2-!^) + /iz,(^)^4 



D(x) 1 f dz ^ /' X 



ln(l - z) 



1 - z 



(1-z) 



- (1 + z) In 



zj yii-z)+ 

-2(l + z)ln(l-z) 

-I- 

_2 z 1 

(1 -z)+r2 + z(l-r2) I 



-(l + z)-5(l-z) 



(23) 



The infra-red and collinear divergences explicit in the result for the integrated dipole Eq. (j23| 
must cancel with the virtual corrections and the MS renormalization of the fragmentation function 
The sum of the leading order decay rate and the one-loop virtual correction reads 



dr^°^+dr^''^^|X2|"/virt d$2, 



where 

^virt = D{x) 
The functions C'o,i read [1] 



a^CF(l - r2)-2<^ /47r^' 



27rr(l - e) 



Co 



2 H-2r2 



(24) 
(25) 



C„^-l-|-li±^4-2Li.(.V21n(l-.^) 



Taking the sum of /dip and /virt , we obtain 



ln(l - r^) 



(26) 



/dip + /virt ^ D{x) + 

1 /■ dz 



-D 

z \z 



(f)p„(=) 



27rr(l-e) V 
dz 



-D 

z \ z 



(1-z) 



- (1 + z) In 



+ z(l - r^) 



(1-rV 

"ln(l - z)' 

1- z _ 

2 



/:»(x)P"(r) 
-2(H-z) ln(l -z) 



(1 - z)+ r2 + z(l - r2) 



(27) 
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where Pqq = 2/(1 - z)+ - (1 + z) + 3/2 6{l - z) and, for = 1, 



V{t) 



-^-^-2Li2(.^)-^ln(l-r^) 
2 3 ^ ^ l + 2r2 ^ ^ 



(28) 



We observe that Eq. (1271) contains colhnear divergences. To remove them, the fragmentation 
function D{x) needs to be renormahzed. By convention, we use the MS scheme. We obtain 



/virt + /dip = i?(/i, x) + ^iM^ ( D{p., x)V{r) - In , ^ 



27r 
ln(l - z) 



mj (1 — r^) 



dz 



(1-z) 



(1 + z) In 



1 - z 

+ z(l - P) 



-2(l + z)hi(l-z) 
2 



(1 - z)+ r2 + z(l - r2) 



(29) 



Equation (I29p contains everything that is needed to compute the contribution of the virtual 
corrections and the integrated dipoles to the decay rate t — > l^v + B + X. These results should 
be supplemented with the contribution of the real emission matrix elements, described by Eq. 
Combining Eq. (I^Hl) and Eq. we can compute ©(as) correction to the fully differential rate 
for t — > l^iy + B + X. We then interface the corrections to the decay, that we just described, 
with the production process, in the spirit of Ref. This allows us to get a description of 

— >■ (t — 7> + & — >■ W+ + J/')/') + (f — >■ W~ + b) at leading and next-to-leading order, including 
the possibility to apply kinematic cuts to the final state particles. 

2.2. The fragmentation function 

The NLO QCD calculation described in the previous Section leads to radiative corrections 
enhanced by the logarithm of the ratio of the top quark mass and the factorization scale /i. 
We can choose mt to get rid of the logarithmically enhanced terms in the short-distance 
partonic decay rate (cf. Eq. ([29])). However, by doing that, we face the challenge of evaluating the 
fragmentation function D{fj,,x) at a high value of the factorization scale in spite of the fact that 
b B fragmentation is, intrinsically, the low-scale phenomenon. 

The standard way to deal with the problem is to use the Altarelli-Parisi (AP) equation 



>dDh 



E 



dz 



-Pj 



\ z 



D 



(z,/x). 



(30) 



to evolve the fragmentation function to the required values of the factorization scale fi ^ nit. 
For the purpose of the NLO calculation, we include 0{as) and 0{a1) contributions to the AP 
evolution kernel which leads to a resummation of the leading and next-to-leading logarithms of 
the ratio of the factorization scale and the fe-quark mass. Similar to what was done in the previous 
studies, we neglect all off-diagonal contributions to the evolution equation Eq. (l30l) and only keep 
there terms proportional to Ptb splitting function. 

Solution of the AP equation requires an initial condition, which is to say that D^^b needs 
to be known for some value of the factorization scale /^o- Traditionally, this is accomplished by 



fitting the fragmentation function at the scale /ip ^ "nib to data on e+e~ — bb 25|, |26[. Since 
/LtQ mi, is a perturbative scale, we may attempt to completely factorize perturbative and non- 
perturbative contributions, by writing the heavy quark fragmentation function as a convolution 
of the perturbative fragmentation function Db{fi,x) and the "non-perturbative" fragmentation 
function I?np(a;) [21 1 



b^B 



= J ^Db{^,,ODn, (1 



(31) 
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The perturbative fragmentation function receives contributions from momenta comparable to the 
6- quark mass and is therefore computable in perturbation theory. At NLO QCD, the result reads 



21| 



Dbifi, X) = 6(1 - x) H 

Zn 



Lt^ln(^]-2\og{l-x)-l +0{ai). (32) 



+ 



Note that the expansion parameter in Eq. p2p is as log /i/mfo. This observation makes it clear 
that D},(fj,, x) has to be evaluated at the scale /i ^ mi,. 

In Ref. [2^ D^pix) was determined from fits to the e+e^ — > bb data by ALEPH and SLD 



collaborations [25|, |26|. For numerical calculations, we use results reported in that reference. We 



restrict our attention to two types of non-perturbative fragmentation functions 

n _ ( x'^il~xr/B{a + l,(3 + l), 
"P~\ {1 + S){2 + S)il ~ x)x^ . ^'^■^> 



It was shown in Ref. [221 that the following values of the parameters 



a = 0.66 ±0.13, /3 = 12.39 ± 1.04, (5 = 14.97 ±0.44, (34) 

lead to a good fit to the ALEPH data provided that no soft gluon resummation is applied to 
the fragmentation functior0. We will use the range of parameters shown above to estimate the 
sensitivity of the extracted value of the top quark mass to the employed model of the heavy quark 
fragmentation function. We solve the AP evolution equation in a standard way by applying the 
Mellin transform since the AP equation becomes ordinary differential equation in the Mellin space. 
The results that are required to perform the Mellin transform can be found in Refs. (21I. [23|. 

2.3. Results: niBi distribution in top quark decays 

In this Section, we discuss our results for the rriBi spectrum, as obtained within the perturbative 
fragmentation function framework. We ignore all the subtleties associated with the heavy quark 
production mechanism and study the invariant mass of the lepton and the _B-meson as produced 
by the top quark decaying in isolation. We quote results at leading and next-to-leading order, 
but we need to clarify what we mean by that. Indeed, a fragmentation function, is extracted 
from data on e+e~ annihilation to _B-hadrons, using a short-distance function for e^e~ — 66, 
computed through a particular order in perturbative QCD. Therefore, if we change the short- 
distance function by truncating it to leading order, we get a different fragmentation function. 
This phenomenon is well-known from studies of parton distribution functions that do change 
from one order in perturbation theory to the other. Unfortunately, information on how non- 
perturbative fragmentation functions change when perturbative predictions for e~'"e~ — )■ 66 are 
truncated at leading order are not available to us, so that for our leading order calculation we 
use the same non-perturbative fragmentation function D^p, Eq. (j33[) . as in next-to- leading order 
computation. However, for leading order computations, we neglect all the 0{as) corrections to 
partonic decay rate of the top quark and the initial condition Db, Eq. (I32|) . and we solve the 
AP evolution equation and compute the evolution of the strong coupling constant in the leading 
logarithmic approximation. To obtain numerical results reported below, we use as{Mz) — 0.130 
and as{Mz) = 0.118 for leading and next-to-leading computations, respectively. 

In Tables 1 and 2, we show average values of the invariant mass of the _B-meson and the 
lepton {msi) and the dispersion cTms, of the msi distribution at leading and next-to-leading 
order in perturbative QCD. To arrive at those results, we calculate {msi) for three values of the 
renormalization and factorization scales mt,mt/2,mt/4:, changing them independently. We also 
use two different fragmentation functions, as explained in the previous Section. For each parameter 
that one needs to describe the fragmentation function, we do a calculation for its central value 



* If soft gluon resummation in the perturbative fragmentation function is employed, the preferred values of a, p 
and 5 change, see Ref. |22| | for details. 
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nit 


{iriBi), LO 


{niBi), NLO 




(7„,^,, NLO 


171 


73.51 ±1.87 


76.03 ±0.61 


31.46 ±0.12 


29.21 ±0.29 


173 


74.71 ±1.90 


77.24 ±0.62 


31.92 ±0.12 


29.63 ±0.29 


175 


75.91 ±1.93 


78.44 ±0.63 


32.37±0.13 


30.04 ±0.30 


177 


77.10 ± 1.95 


79.64 ±0.63 


32.82 ±0.13 


30.46 ±0.30 


179 


78.29 ± 1.98 


80.84 ±0.64 


33.26 ±0.13 


30.87 ±0.30 



Table 1 : The estimate of the average value of the B-meson-lepton invariant mass and its dispersion at leading and 
next-to-leading order, in dependence of the top quark mass. The top quark mass and all the results are in GeV. 
Decay of an isolated top quark is considered. 



nit 


(m_B/), LO 


{niBi), NLO 


'^msn LO 


(JmBU NLO 


171 


87.51 ±1.04 


87.20 ±0.43 


22.17±0.23 


21.28 ±0.17 


173 


88.53 ± 1.07 


88.22 ±0.43 


22.68 ±0.24 


21.77 ±0.17 


175 


89.56 ± 1.10 


89.25 ±0.44 


23.19 ±0.24 


22.25 ±0.18 


177 


90.58 ± 1.13 


90.29 ±0.45 


23.69 ±0.24 


22.73 ±0.18 


179 


91.61 ± 1.15 


91.32 ±0.46 


24.20 ±0.24 


23.22 ±0.18 



Table 2: The average value of the invariant mass {mgi) and its dispersion, evaluated with the cut on the invariant 
mass mBi > 50 GeV. The top quark masses and all the results are in GeV. Decay of an isolated top quark is 
considered. 



and for the central value shifted by plus/minus the error quoted for that parameter. As the result, 
we obtain 108 values of {mriBi) and CTniB, for each of the input values of the top quark mass. We 
calculate the mean and the error from these samples of 108 numbers for both (tob;) and <7mBi 
at leading and next-to-leading order. In Tables 1 and 2 those results are shown; the difference 
between the two Tables is that an additional constraint ttibi > 50 GeV is employed to obtain 
results in Table 2. 

There are two immediate comments that one can make about those results. First, we observe 
that NLO QCD corrections to (mBi) strongly depend on the applied cut on the invariant mass. 
For example, if no such cut is applied, the shift from leading to next-to-leading order in (mBi) 
is about 2.5 GeV, whereas if a 50 GeV cut is applied, {niBi) shifts by —0.3 GeV. Second, the 
uncertainty in {tjibi) decreases by a factor between two and three, when NLO QCD effects are 
included, indicating their importance for the high-precision top quark mass measurement. 

It is interesting to compare the results of the computation reported in this paper with the 
previous analysis where parton shower event generators were employed 18. [l^. We note that 



results of those two references are not consistent; the reason is explained in Ref. [l9|. We will 
therefore compare to the results in Ref. ^1^, where (niBi) and higher moments of Bl invariant 
mass distribution are computed using HERWIG and PYTHIA. The B ^ b fragmentation functions 
were fitted in Ref. [lij to reproduce i3-meson energy spectra in e+e^ annihilation. Systematic 
0(1 GeV) differences in values of (tob;) obtained with PYTHIA and HERWIG were observed in 
[l9j, with PYTHIA results being lower. 

We find that the NLO QCD result for {niBi) and (Tmg, are close to the results obtained with 
parton showers. Nevertheless, the difference is not negligible, given the expected precision of the 
top quark mass measurement. By comparing our results with that of Ref. [l9| . we find that the 
average values of {ruBi) computed through NLO QCD is about 2.4 GeV lower than {niBi) obtained 
with HERWIG and only 1.1 GeV lower than {niBi) obtained PYTHIA. On the other hand, the 
dispersion aBi that we compute through NLO QCD, differs by 1 GeV from PYTHIA and by 
0.5 GeV from HERWIG results. 

The results for {ttibi) displayed in Tables 1 and 2 can be described by a linear function of the 
top quark mass. We present the results of such a fit for the two cases - with and without a cut 
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m, I GeV ] m, [ GeV ] 

Figure 1: Results of the linear fit to (mgi)^^'-^ are shown. Left panel - no cut on nigi is applied. Right panel - 
'tis; > 50 GeV cut is applied. In both cases, decays of isolated top quarks are considered. 



on niBi in Fig- [2 We find 

= O.eOlmt - 26.7 GeV, = 0.004; (35) 

{mBi)^t^%o GeV - 0.516mt - 0.96 GeV, (5™^ = 0.006, (36) 

where Si-ms is the root mean square (rms) of the residuals of the linear fit. It is clear from the 
value of (Jims that the linear fit works very well. It is straightforward to translate the results of the 
linear fit shown in Eqs. P5I36I) to an estimate of the error on the top quark mass. Indeed, suppose 
that a typical uncertainty of the measured value of {msi) is 0.4 GeV. The slopes in Eqs. (|35l36p 
then imply that the corresponding error in the top quark mass mt is about 0.8 GeV. On the other 
hand, assuming perfect measurement of (mBi), we find that theoretical uncertainties in (mBi) 
shown in Table 1 and the value of the slope of the linear fit lead to a 1 GeV uncertainty in the 
extracted value of the top quark mass. The errors on the top quark mass that follow from the 
NLO QCD computation are similar to differences between PYTHIA, HERWIG and NLO QCD. 
To show this, we quote results of a fit to {niBi) obtained with PYTHIA and HERWIG in Rcf. [l^ 

(msi)pythia = 0.59 mt - 24.11 GeV, (mi3i)Hcrwig = 0.61 rut - 25.31 GeV. (37) 

It is clear from the comparison of the fits Eqs. p7l35p that NLO QCD results and parton shower 
results are close but not identical and these differences are essential. Indeed, we note that a slope 
difference between parton showers and NLO QCD is about 0.01. Although such slope difference 
may look insignificant, it leads to 0(3 GeV) shift in the reconstructed value of the top quark. 
Hence, parton showers are insufficient for measurements of the top quark mass with a precision 
higher than a few GeV. On the contrary, it follows from Eqs. P5l36p that NLO QCD computations 
lead to results with small uncertainties that can be estimated in a systematic way and, perhaps, 
be even further improved. 

2.4- Results: niBi distribution in pp ^ {t ^ + 6 + i?) + (i W~ + b) 

In this Section, we consider production of _B-mesons through the fragmentation of 6-quarks in 
top decays but, in contrast to the previous Section, we include the full production process through 
next-to-leading order in perturbative QCD. To claim that the NLO QCD computation can do a 
good job in describing (mBi) in reality, it is very important to have full production and decay 
chain included. 

To this end, we consider top quark pair production in pp collisions at ^/s = 14 TeV and focus 
on the lepton -I- jets decay channel. We note that NLO QCD corrections to the decay W — >■ qq' 
need to be interfaced with pp — >■ tt production process, to describe lepton -I- jets channel through 
NLO QCD. We require that there are at least four jets in the event. We include the S-meson 
and the non-S-meson remnant of the fragmenting 6-quark into a list of proto-jets that are passed 
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nit 


{iriBi), LO 


{niBi), NLO 




(7„,^,, NLO 


171 


77.07 ± 1.92 


76.75 ± 1.12 


30.60 ±0.13 


28.41 ±0.36 


173 


78.34 ± 1.93 


77.92 ± 1.09 


31.01 ±0.14 


28.72 ±0.31 


175 


79.63 ± 1.98 


79.31 ±1.04 


31.46 ±0.14 


29.12 ±0.18 


177 


80.91 ±2.03 


80.55 ± 1.05 


31.83 ±0.15 


29.48 ±0.13 


179 


82.16 ±2.04 


81.80 ± 1.04 


32.24±0.16 


29.83 ±0.13 



Table 3: The average values of the invariant mass (msi) and the dispersion in case where all the cuts on the final 
state particles are applied. The top quark masses and all the results are in GeV. See text for details. 



to the jet reconstruction algorithm. We employ k± jet algorithm with R = 0.5 and the four- 
momentum recombination scheme. All reconstructed jets and the positron from the decay 
are required to have transverse momenta in excess of 20 GeMj- The scalar sum of the transverse 
momenta of all jets in the event should exceed 100 GeV [20|. For the sake of simplicity, we do 
not impose any other kinematic constraints including cuts on the missing energy and the lepton 
isolation cuts. In addition, we do not consider combinatorial backgrounds, assuming that the 
correct pairingbetween a lepton and a _B-meson can be established. Finally, similar to what was 
done in Ref. [l3|, throughout this paper we consider intermediate to p q uarks to be on the mass- 
shell and we do not include the so-called non-factorizable corrections [27| . For observables that we 
study in this paper, this is a good approximation since we, effectively, integrate over the invariant 
masses of each of the top quarks. 

In contrast to the previous Section, we do not change the parameters of the fragmentation 
functions, fixing them to their central values, see Eq. (IMl) . We use CTEQ parton distribution 
functions [2^ [2^ in the analysis. For each input value of we compute (mBi) for three values 
of the renormalization and (pdf)-factorization scales = fip — ['Ti(/4, mt/2, rrit] and for three 
values of the factorization scale in the b ^ B fragmentation function [mt/4, mt/2, mt] and for two 
different types of fragmentation functions Eq. (p4| . As the result, for each value of rrit we have 
eighteen values of {niBi) and (JmBi- Although this is not an extensive scan of the parameter space, 
it gives a sense of theoretical uncertainties in {vriBi) provided that realistic production mechanism 
is employecll. Computing the mean and the error, we arrive at the results shown in Table 3. We 
see that the NLO QCD effects in this case are quite small and negative, which is similar to the 
case oiruBi > 50 GeV cut discussed earlier. It is also clear that the uncertainty in tubi decreases 
when NLO QCD corrections are included. By comparing results in Table 1 and in Table 3, we see 
that effects of kinematic cuts on {niBi) are more important at leading order, where the average 
value of niBi can shift by 3.5 GeV. On the other hand, at NLO, a typical shift is of the order of 
0.8 GeV and, therefore, is less dramatic. Comparison of Tables 1 and 3 shows that the uncertainty 
of (rriBi) at NLO nearly doubles if production mechanism is taken into account. Performing the 
fit, we obtain (see Figl2|) 

{mBi)ptd = 0-6365 mt - 32.12 GeV, <5„ns = 0.053. (38) 

Comparing this result with Eq. (|35|) . we find a significant change in both the slope and the 
constant part. This demonstrates that (mBi) depends in a non-trivial way on the production 
mechanism, because of kinematic cuts applied to top quark decay products and additional jets in 
the production process. It follows from Eq. (|38|) and the uncertainties of {niBi) shown in Table 3 
that the theoretical error on the extracted value of mt is close to 1.5 GeV. 



standard argument [T^. Il9ll that {vagi) involves a Lorentz invariant product of the two four-vector and, 
therefore, does not depend on the production mechanism is not applicable once cuts on the transverse momenta are 
applied. Such cuts are only invariant under restricted class of Lorentz transformations - boosts along the collision 
axis. 

®We have checked that if we only change parameters that are related to the decay process, we find the 0(0.5 GeV) 
uncertainty in (mgi), similar to Tables 1,2. 
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Figure 2: Result of the linear fit to (mgi)^^^ is shown, with all kinematic cuts on the final state particles applied. 
See text for details. 



3. Dilepton channel 

In the previous Section we saw that top quark decays to final states with identified hadrons 
provide an interesting way to determine the top quark mass. In this Section we study inclusive 
final states. We focus on the case where the top and the anti-top quarks decay semileptonically, 
e.g. t W^b — l^vb. We study the kinematic distribution of an invariant mass of a 6-jet 
and a lepton, and the distributions of the sum of energies of the two leptons and the two 6-jets. 
We employ the NLO QCD corrections to top quark pair production and decay, as computed in 
Ref. [13j. Throughout this Section, the center-of-mass energy of proton collisions is 14 TeV. 

We begin by summarizing the kinematic cuts that are employed to identify dilepton events 
Leptons are required to be central [yy'l < 2.5 and have large transverse momentum > 25 GeV. 
There should be missing energy in the event, Ef^^^ > 40 GeV. The jet transverse momentum cut 
is p±j > 25 GeV. We employ the k± jet algorithm with R = 0.4. 

3.1. Invariant mass of a lepton and a b-jet 

It is pointed out in Ref. Q that an average value of the invariant mass squared of a 6-jet and a 
lepton and an average value of the the angle between the lepton and the 6-jet in the W boson 
rest frame, can be used to construct an estimator of the top quark mass. The estimator reads 

ML -rnl.+ . (39) 

1 - [cos ff lb) 

To see that this is a good estimator, we note that for the top quark decay computed at leading 
order in perturbative QCD and without any restrictions on the final state Mcst equals to mt 

{ml) = ^^^l^^{l-{coseib)), (coseib) ^ ^ML = ml (40) 



w 



In reality Meat is not equal to rrit for a variety of reasons including i) kinematic cuts required to 
identify the dilepton events; ii) effects of higher order QCD corrections; iii) impossibility to choose 
the "correct" pair of a lepton and a 6-jet and iv) the experimental issues with 6-jet misidentification 



and the jet energy resolution. The computation reported in Ref. [13[ allows us to calculate 
within the framework of perturbative QCD, accounting for the points i)-iii) exactly. 

We point out that the computation of NLO QCD corrections to pp — >■ tt process reported in 
[l3| includes exact spin correlations, one-loop effects in top quark decays and allows arbitrary 
constraints to be imposed on top quark decay products. These features are crucial for reproducing 
experimental procedures. Indeed, experimentally, it is not possible to determine the charge of the 
6-jet. Hence, it is unclear which of the two 6-jets should be combined with the chosen, definite- 
sign, lepton. For the purpose of mib reconstruction, one pairs the lepton with the 6-jet that gives 
the smallest mib value [J]. The parameter {cos 9ib) in Eq. ([39]) is not measured and should be 
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Figure 3: The invariant mass distribution of the lepton and the 6-jet. Note that the lepton and the 6-jet do not 
necessarily come from the decay of the same top quark, see text. The left panel shows the scale uncertainty bands 
for /iji = fip = [0.5mi, 0.75r?it, mt, 1.25mt]. The right panel shows two NLO normalized mn, distributions for 
mt = 171 GeV and mt = 179 GeV. 



estimated theoretically. We have also chosen to calculate {cos On,) for the 6-jet that minimizes 
the invariant mass mu, since in this case, there is a partial compensation of incorrect assignments 
between the numerator and the denominator in Eq. (jJH). As the result, Mest becomes closer to 
the input value mt as compared to the case when "correct" pairing of the 5-jet and the lepton 
is chosen to calculate {cosBw) in Eq. (|39|) . It is argued in Ref. that with 10 fb"""^ integrated 
luminosity, the statistical and systematic uncertainties in the top quark mass of about 1 GeV each 
can be achieved from (mf^) measurement. 

To assess how realistic those uncertainties are, we consider five different values of the top quark 
mass mt — [171, 173, 175, 177, 179] GeV. For each of these m* values, we compute Afcst for four 
values of the renormalization and the factorization scales Mr =_Mf — [O.Stoj, 0.75mt, mj, 1.25mt] 
and for two sets of parton distribution functions CTEQ [3, [l^ and MRST [s^l- We use the 
mean value and the standard deviation of these eight values to compute central value of Most and 
its error. Clearly, by no means this is an exhaustive scan through the parameter spacc[ll but it 
gives us an idea of the uncertainties on the theoretical side. Examples of mu, distributions and 
the results of the calculation are shown in Figs. I3I4I The uncertainties on Mest do not depend on 
mt in significant way; they are 0.1 (0.2) GeV at leading and next-to-leading order, respectively. 
Performing the linear fit, we find 

Mg^o = 0.8262TOt + 23.22 GeV, M^^J^o = 0.7850mt + 28.70 GeV. (41) 

The quality of the linear fit is very good; for example, the root mean square of the residuals of 
the NLO fit is 5rms = 0.032. It is instructive that the analysis of this observable at leading order 
shows stronger correlation between m* and Mest than at next-to-leading order. In addition, the 
theoretical uncertainty in Mest increases when NLO QCD corrections are included. The primary 
reason for the increased uncertainty is stronger dependence of Mest on the renormalization and 
factorization scales at NLO. This feature can be understood by considering the situation where 
no phase-space cuts are applied and where all the assignments of a lepton and a &-jet are done 
correctly. In this case, as follows from the discussion at the beginning of this Section, the estimator 
equals to the top quark mass regardless of the renormalization and factorization scales and the 
chosen parton distribution functions. At next-to-leading order, this is not true anymore because 
of the gluon radiation in top decay that is sensitive to the value of the strong coupling constant 
and, hence, to the renormalization scale. We note that we observe a very weak dependence of Mest 
on parton distribution functions which implies that even with the phase-space cuts and incorrect 
pairing, this variable is primarily sensitive to top quark decays rather than to top quark production 



^For example, one can and perhaps should use different renormalization scales to compute numerator and 
denominator in Eq. I I39II . to get a better idea of the scale uncertainties in Mest- 
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mechanism. 

Finally, we can use Eq. (j4T|) and Fig. |4] to estimate uncertainty in mt that can be achieved by 
measuring Mest with infinite precision. Since, as follows from Fig. SI the uncertainty in Most is 
0.2 GeV at NLO and given the slope of about 0.8 in Eq. (|4T|). we find the minimal uncertainty in 
the extracted value of rrit to be close to 0.25 GeV. We note that this result does not include the 
6-quark fragmentation uncertainty and the jet scale uncertainty, estimated to be 0.7 and 0.6 GeV, 
respectively, in Ref. 0]. 

3.2. Sum of energies of the two leptons from top quark decays 

Another observable that we consider is the sum of the energies of the two leptons, Ei-^ + Ei^, 
in the laboratory frame. Lepton energies in the laboratory frame can be easily measured and they 
are free from jet energy scale uncertainties that are important sources of errors, if the top quark 
mass is reconstructed from hadronic final states. The important question is whether or not the 
average valu(^ of the sum of lepton energies is correlated with the top quark mass at the parton 
level and how well such correlation can be described by perturbative QCD. 

The corresponding distributions are shown in Fig. [5l There we display Ei-^^ + Ei^ computed 
through leading and next-to-leading order in perturbative QCD for mt = 175 GeV, as well the 
NLO QCD distributions in Ei^ -f Ei^ for rrit = 171 GeV and mt = 179 GeV. To compute the 
mean value of {Ei-^ + i?;^.) we consider the same range of the renormalization and factorization 
scales and the two sets of parton distribution functions, as in the previous Section. The results of 
the calculation are shown in Fig. [6l Performing a linear fit, we find 

{Ei,+Ei,)w = 0.645mt + 120.6 GeV, <5„ns = 0.08; 

{El, + ^,,)nlo = 0.670mt + 114.4 GeV, 5,,^, = 0.07. (42) 

The results of the linear fit are displayed in Fig. [51 Theoretical errors on {Ei, +E12) are independent 
of the top mass; they are 1.7 GeV at leading order and 1 GeV at next-to-leading order. Combining 
information about the slope in Eq. (|42|) with the theoretical uncertainty on {Ei-, +E12), we conclude 
that the ultimate uncertainty in mt that can be achieved by studying this observable is close to 
1.5 GeV. It is interesting to point out that, in this case, both the scale dependence of the NLO 
result and the diff'erence between MRST and CTEQ parton distribution functions are the two 
important sources of the uncertainty. 

3.3. Sum of jet energies 

Another observable that was discussed [3| in connection with the top quark mass measurement 
is the sum of energies of the two hardest jets in the laboratory frame. Similar to the lepton energies 



*One can ask the same question about the shape of the distribution but such discussion is outside the scope of 
this paper. 



14 




[GeV] [GeV] 

Figure 5; Left panel: normalized distribution of the sum of lepton energies at leading and next-to-leading order 
calculated for mt = 175 GeV. The renormalization and factorization scales are set to mt and the MRST (left panel) 
and CTEQ (right panel) parton distribution functions set is used. Note a shift in the position of the maximum of 
this distribution. Right panel: normalized distributions of the sum of lepton energies at next-to-leading order, for 
mt = 171 GeV and mt = 179 GeV. 
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Figure 6: Results of a linear fit to the sum of the average energy of the two leptons at leading and next-to-leading 
order. 



just discussed, the shape of the distribution is an observable that is to be fitted; this is beyond the 
scope of the present paper. Here, we limit ourselves to the discussion of average values. Instead of 
considering the two hardest jets in the event, we found it more useful to take the sum of energies 
of the two 6-jets. 

We consider the distribution of the sum of energies of the two 6-jets. We expect that this 
distribution is strongly correlated with the top quark mass, since 6-quarks originate directly from 
top decays. The results of the calculation are shown in Figs.[71 Performing a linear fit, we obtain 
the correlation between the average value of the two 6-jets and the top quark mass 

{Ebji + Ebj2ho = 2.18TOt - 42.2 GeV, (5™s = 0.02; 

{Ebji + Ebj2)NLO = 2.09mt - 29.2 GeV, (5„„s = 0.05. (43) 

The results of the linear fit together with theoretical uncertainties in {Ebji -f Ebj2) are shown 
in Fig. m These uncertainties are 2.6 GeV at LO and 2.4 GeV at NLO; they do not exhibit a 
strong dependence on the top quark mass. Interestingly, inclusion of NLO QCD corrections makes 
the correlation between {Ebji + i?6j2)LO and mt weaker. However, the correlation is still quite 
strong. If we assume that energies of 6-jets can be measured infinitely accurately, the irreducible 
uncertainty on the determination of the top quark mass from {Ebji +Ebj2) becomes only 1.2 GeV. 
Of course, the main issue here is to understand how well 6-jet energies can actually be measured; 
this issue will be at the center of the experimental studies at the LHC. 
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Figure 7: Normalized distribution of the sum of energies of two b-jets -Bbjet = ^bjl + ^bj2 next-to-leading order, 
calculated for mt = 171 GeV and mt = 179 GeV. We set renormalization and factorization scales to mt and use 
CTEQ parton distribution functions. 
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Figure 8: Results of a linear fit to the sum of the average energy of the two 6-jets at leading and next-to-leading 
order. 



4. Conclusion 

Determination of the top quark mass with high precision is an important part of the top 
quark physics program at the LHC. It is expected, that a variety of methods will be employed by 
ATLAS and CMS collaborations to measure the top quark mass. Some of those methods involve 
studies of the correlations between the value of the top quark mass and the kinematics of the 
top quark decay products. As the results of those studies, it is often claimed that the top quark 
mass can be determined with 0(1%) uncertainty at the LHC, but it is not clear whether or not 
these uncertainty estimates can be trusted. Indeed, all such studies employ PYTHIA or HERWIG 
parton shower event generators to describe top quark pair production and decay; however, no 
parton shower is designed to handle this level of precision. An interesting question therefore is 
to estimate, in a parton-shower-independent way, the uncertainty on mt that can be achieved in 
forthcoming LHC experiments. 

The goal of this paper is to address this question by computing correlations between the top 
quark mass and the values of some kinematic observables through NLO QCD. Among other things, 
we consider NLO QCD corrections to the invariant mass of a i?-meson and a lepton from top quark 
decays which is considered to be one of the most accurate ways to determine the top quark mass. 
Such computation is rather unusual in the context of NLO QCD calculations since it refers to the 
final state with an identified hadron. 

In general, we find that parton shower event generators do a good job in estimating both the 
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central value and the uncertainty in the top quark mass that can be achieved. However, as can 
be seen from the discussion of the average value of the invariant mass of the i?-meson and the 
lepton, NLO QCD computations give both, a more accurate central values and an estimate of 
the uncertainty that can be trusted. Both of these features are important if we want to use the 
measured value of the top quark mass with confidence, to constrain physics beyond the Standard 
Model through precision measurements. 
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